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The melting transition in the hard-disk system is considered. Non-equilibrium relaxation 
analysis of the six-fold bond-orientational order parameter has been carried out. The critical 
(^ ' point between the hexatic and the fluid phase is determined on the basis of the dynamic 

fT^ , scaling hypothesis. The value of the critical exponent 77 is determined from the fluctuation of 

the order parameter at the criticality as 77 = 0.25(2) which is consistent with the prediction 
by the Kosterlitz-Thouless-Halperin-Nelson- Young theory. 
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In 1957, Alder et al. reported the melting transition in the hard-sphere sys- 
tem. ^^ The fact that particles only with repulsive interaction involve the phase 
transition has surprised many researchers. This transition was also confirmed in the 
Cd ' hard-disk system,^-' and these kinds of transitions are now often called the Alder 

P^ ■ transition. The Alder transition seems to be the usual first-order type since it is 

F^ i characterized by a van der Waals loop like behavior in the pressure-density diagram 

^ ' (see Fig. 1). However, Mermin has proved that the long-range correlation of the 

O . translational order is inhibited, ^^ and therefore, the two-dimensional melting cannot 

be the usual order-disorder phase transition. Halperin, Nelson and Young proposed 
the theory of the two-dimensional melting based on the unbinding mechanism of 
^ . the defects.^'' ^^ They have suggested that the two-dimensional melting involves 

t — I the Kosterlitz-Thouless transition twice. This theory is now often referred to the 






Kosterlitz-Thouless-Halperin-Nelson- Young (KTHNY) theory. Competing theory 

predicting the first-order transition was proposed by Chui.®-* Many experimental 

and numerical studies have followed in order to clarify the nature of the transition. 

(^ ' In spite of these efforts, the strong finite-size effect has kept us from fixing this 

00 , problem and the nature of the two-dimensional melting remains unanswered.^) 

In this paper, an overview of the two-dimensional melting is given first. After 
the Mermin's theorem and the KTHNY theory are described, a new method, which 
^ . utilizes the non-equilibrium relaxation behavior of the order parameter, is introduced 

H I in order to avoid the finite-size effect. Conclusions and some recent developments 

are summarized in the final section. 
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§2. The KTHNY theory 



2.1. Translational order 



Roughly speaking, states of matter can be classified into two phases, fluid and 
solid. Fluid is a substance which deforms under an infinitesimal shear stress, while 
solid requires finite stress. In this definition, gas and liquid are not discriminated 
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Fig. 1. The density-pressure diagram of the particle-system only with the repulsive interaction. 
The melting behavior can be confirmed because of the existence of the van der Waals loop. 
Note that, this figure is exaggerated. The domain of the loop becomes very small as the system 
size becomes larger. 



and the both are classified into fluid. Crystal is a solid with the translational order. 
Note that, the above definitions are not absolute, since there are many substances 
which are difficult to be classified such as glasses and liquid crystals, etc. Instead, it 
is sometimes appropriate to define that the Alder transition itself distinguishes fluid 
from solid. In this section, we review briefly the deflnition of the translational order 
in the crystal, and how it is inhibited by the Mermin's theorem. 

Consider a system with N particles. The density distribution function p{r) can 
be defined as 



P(r) 



N 

E 



(5(r 



with the positions of the particles r^ . The fourier transform of the density p^ is given 
by 



Pk = / dr exp (ik • r)p(r) 

N 
= ^exp(ik-rj). 



f2-l) 



(2-2) 



If the system has the translational order, at least one reciprocal vector K. With 
this wave vector, the value of Eq. (2-2) becomes 0{N). Therefore, the translational 
order parameter of this system is defined by 



^ 



1 ^ 
— ^exp(iK-rj). 



(2-3) 



Mermin has proved that the translational order parameter •P' cannot exhibit long 
range correlation in the two-dimensional system, provided that the the two-body 
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potential <P{r) satisfies the following conditions 

(p{r) < r^^+l^l (r -^ oo), 



^(r) > r-2+kl (r ^ 0). ^^"^^ 

The translational order 'F of two-dimensional system decreases faster than logarith- 
mically as the system size increases, and will vanish in the thermodynamic limit. 
There are no such constraints in the three-dimensional system, therefore, the nature 
of the Alder transition in two-dimensional system is fundamentally different from 
that of the three-dimensional system. 

The following points are worth noting. The Mermin's theorem cannot be applied 
to the self-gravity system, since the system does not satisfy the condition (2-4). It 
is not trivial whether the theorem can be applied to the hard-particles, since the 
hard-core potential is not differentiable. 

2.2. Bond-orientational order 

The constraint on the translational order in the two-dimensional system is so 
week that some kinds of orderings are allowed to exist. One of such the orderings is 
the bond-orientational order which plays an important role in the KTHNY theory. 
The six-fold bond-orientational order (J)q is defined as 

(/>6 = (exp(6i^)) (2-5) 

with the angle 9 between a fixed axis and the bond connecting neighboring particles. 
The average is taken for all pairs of neighboring particles. The parameter (pQ becomes 
1 when all particles are located on the points of the hexagonal grid, and it becomes 
when the particle location is completely disordered. Therefore (J)q describes how 
close the system is to the perfect hexagonal packing. The neighbors in an off-lattice 
model are defined with the Voronoi construction. 

Unlike the translational order parameter, the bond-orientational order parameter 
does not change its value under uniform dilation of the system. The long-range 
correlation of the bond-orientational order is not inhibited in the two-dimensional 
system while the translational order is. 

2.3. Two kinds of defects 

The two-dimensional classical XY spin model cannot have the spontaneous mag- 
netization because of the Mermin-Wagner's theorem. The system, however, involves 
the phase transition involving divergence of the susceptibility.^-* This phase tran- 
sition is called Kosterlitz-Thouless (KT) transition which has some characteristic 
properties such as the low temperature phase with power-low correlation and the es- 
sential singularity at the critical point. The mechanism of this transition is explained 
by the unbinding of oppositely charged topological defects. 

Similar unbinding mechanism can be applied to the two-dimensional melting. 
In the two-dimensional crystal, there are two kinds of defects which are dislocation 
and disclination. Schematic drawings of the defects are shown in Fig. 2. While 
the dislocation destroys the translational ordering, it does not destroy the bond- 
orientational ordering. The disclination destroys the both ordering, and therefore, its 
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influence is global. Halperin, Nelson and Young have explained the two-dimensional 
melting on the basis of the unbinding mechanism of the two kinds of defects and this 
mechanism is now called the KTHNY theory. The KTHNY theory predicts that the 
two kinds of orders involve the KT transition independently. As density decreases, 
solid melts into anisotropic fluid at density pm and the anisotropic fluid becomes 
isotropic at the density pi. Therefore, the solid becomes fluid via new phase which 
is called the hexatic phase. The KTHNY theory is summarized in Table I. 





Fig. 2. Defects in the triangular lattice, (left) A disclination. The system has five-fold symmetry 
because of the disclination while it has normally six-fold. The charge of the disclination is 
defined by subtracting the number of the neighbors at the defect from the average number of 
neighbors, (right) A dislocation. The charge of dislocations is defined by a Burgers vector which 
is a path around the defect fails to close. Note that the dislocation can be described by the two 
disclination with opposite charges. In this figure, the dislocation consists of the positive (open 
circle) and negative (solid circle) disclinations. 



Phase 


Solid 


Hexatic 


Fluid 


Density 


P> pm 


Pm > p > Pi 


Pi> P 


Dislocation 


Pair 


free 


Free 


Disclination 


Quartet 


Pair 


Free 


translational order 


quasi- long-range 


short-range 


short-range 


Bond-orientational order 


long-range order 


quasi- long-range 


short-range 



Table L Summary of the KTHNY theory. There is the hexatic phase between the solid and the 
fiuid phases. The quasi-long-range order means that the order decays to zero as a power of the 
distance. The theory predicts that two order parameters involve the KT transition indepen- 
dently. 



2.4. First order vs. Continuous transition 

The KTHNY theory depends on some assumptions. One of the important as- 
sumptions is that the defects may be excited tenuously and uniformly. If the defects 
concentrate locally, different arguments are required. Chui calculated the free energy 
on the basis of the collective excitation of the defects, and has concluded the first 
order transition.^) Whether defects are excited uniformly or collectively is depends 
on the value of the core-energy of the effective Hamiltonian for defects.^) The core- 
energy plays a role of the chemical potential. The KTHNY theory is justified for the 
large value of the core-energy, and Chui's theory is justified in reverse. ^'^^ The value 
of core-energy depends on the details of the system, and it is difficult to determine its 
value. ^""^^ Additionally, it is not trivial whether we can apply the elastic theory to the 
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hard-disk system while most of theories are based on the elastic theory. Considering 
the simplified Bernal lattice, Kawamura proposed theory without using the elastic 
theory and has concluded the first order transition. ^^^ 

§3. Non-equilibrium relaxation analysis 

3.1. Finite-size effect on relaxation 

Most of the numerical works studying the Alder transition have used the equilib- 
rium Monte Carlo (MC) simulations. While the MC method is popular and has been 
widely used in statistical physics, this method sometimes faces difficulty in achieving 
the equilibrium state for the system with slow relaxation. It is found that the time to 
achieve the equlibrium for the hard-disk system becomes much longer as the system 
size increases. ^^' Therefore, we introduce a new method, called the non-equilibrium 
relaxation (NER) method, in order to avoid the finite-size effect. ^^^ The main idea 
of the NER method is to utilize the relaxation behavior to the equilibrium state. At 
the beginning stage of the relaxation, the behavior of the system is not influenced 
by the finite-size effect. As time advances, the correlation length grows. Therefore, 
the behavior of the system can be regarded as that of the thermodynamic limit until 
the correlation length reaches the system size. In order to see this finite-size effect 
on the relaxation, time evolutions of the order parameter of the system with several 
sizes are shown in Fig. 3. The time evolution of two different systems are equivalent 
for a certain period of time, and this time becomes longer as the size of systems 
becomes larger. Conversely, the relaxation behavior can be regarded as that of the 
system in the thermodynamic limit provided the behavior are shared in two systems 
with different sizes. 
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Fig. 3. Finite-size effect on the relaxation. Time evolutions of the bond-orientational order param- 
eter 06 are shown for the systems of A^ = 1600, 10000 and 14400. (a) Comparison between the 
small and the large systems. The behaviors begin to differ from around t — 2. (b) Comparison 
between larger systems. The behaviors are completely equivalent, and therefore, this behavior 
can be regarded as that in the thermodynamic limit. 
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3.2. Dynamic scaling anaysis 

In order to investigate the phase transition from the relaxation behavior, we 
study the dynamic scahng behavior of the order parameter. We observe the time 
evolution of the bond-orientational order parameter of systems which initial config- 
urations are set to be the perfect hexagonal-packed. The bond-orientational order 
parameter 4>q{p, t) is a function of density and time, and starts relaxing from 1 to the 
value at the equilibrium. Near the criticality, the behavior of 4'Q{p,t) can be scaled 
as 

MP,t) ^ T'^Mt/r), (3-1) 

with relaxation time r and a density-independent parameter A. The scaling function 
is denoted by (pQ. In the KT transition, the correlation length diverges exponentially 
as 

e = 6exp(a/V^) {e = {pi-p)/pi). (3-2) 

According to the dynamics scaling hypothesis, the relaxation time is associated with 
the correlation length as 

r ~ e, (3-3) 

with the dynamic scaling exponent z. Finally, the divergence behavior of the relax- 
ation time r is expected to be 

T{e) = Bexp{A/^/i). (3-4) 

Therefore, we can determine the critical point p; by observing the relaxation of (/)q 
at several densities. ^^^ 

The asymptotic behavior of 4>Q{t) and its fiuctuation at the critical point are 
expected to be 

Mt) ~ i^"/'^ (3-5) 



N 



{Mt)f 



t^/'. (3-6) 



with the critical exponents r] and z, and the dimensionality d = 2}^' Therefore, we 
can determine the values of the critical exponents by observing fiuctuation behavior 
of the order parameter at the criticality. 

3.3. Numerical Results 

We perform the event-driven molecular dynamics (MD) simulation in order to 
observe the time evolutions of the bond-orientational order parameter. The particle 
number N is fixed at 23288 throughout the following results. The periodic boundary 
conditions are taken for both directions of the simulation box. Up to 512 independent 
samples are averaged at each density. 

The time evolutions of the bond-orientational order parameter are shown in 
Fig. 4(a) and the scaling plot is shown in Fig 4(b). From the divergence behavior of 
the relaxation time, the critical point is determined to be pi = 0.893(1). The critical 
exponents at the critical point are determined to be z = 2.5(2) and rj = 0.25(2)."*^^' 
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Fig. 4. (a) Relaxation of the bond-orientational order 06 (i) for various densities from p — 0.878 
to 0.896 Natural logarithms are used for the both axes, (b) Scaling plot of bond-oricntational 
order parameter with appropriately chosen r(e) and A. 




Fig. 5. Comparison of the behavior of the translational order parameter of the systems which time 
evolution are performed by MC and MD. The density is p = 0.890 in (a) and 0.900 in (b). 
Despite many independent samples are averages, the results from MD have strong oscillation. 
The amplitude of the oscillation becomes larger as density increases. 



§4. Summary and Further issues 



We have reviewed the Alder transition and the KTHNY theory. The presented 
numerical results are consistent with the prediction of the KTHNY theory. Espe- 
cially, the obtained value of r/ ~ 1/4 is a strong piece of evidence that the transition 
is of the KT type. We have not studied the transition between the hexatic and the 
solid phase, since the translational order parameter has strong oscillation when the 
time evolution is performed by MD (see Fig. 5). This oscillation is caused by the 
conservation of the momenta, therefore, the results from MC are free from it. It is 
one of the further issues to investigate whether similar technique presented in this 
paper can be applied to MC method. 

The studies on the Alder transition continues today from the both sides of nu- 
merical works and experiments. To the best of our knowledge, the largest simulation 
to date contains 4 x 10^ particles. "^^^ While the scaling analysis of this study has 
supported the KTHNY theory, the possibility of the week first order transition has 
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also suggested from the pressure-density diagram. The experiments of the two- 
dimensional melting have been carried out on various systems, e.g., liquid crystals, 
absorbed gas on graphite, two-dimensional plasma, and so on. Recently, the excita- 
tion behavior of the disclinations have been directly observed in the dusty plasma. ^^' 
In this experiments, the disclinations did not get dissociated even in low enough den- 
sity. Instead, the grain boundary excitation was observed, and therefore, the theory 
by Chui is supported. 

Generally speaking, it is very difficult to distinguish the continuous transition 
from the week first-order transition. Some new methods beyond the analysis of 
the scaling and the pressure-density diagram have been waited such as the level 
spectroscopy method which has achieved success to investigate the KT transition of 
the lattice systems. ^'^' 
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